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Abstract 

We investigate the set a) of positive, trace preserving maps acting on density 
matrices of size A'^, and a sequence of its nested subsets: the sets of maps which are 
b) decomposable, c) completely positive, d) extended by identity impose positive 
partial transpose and e) are superpositive. Working with the Hilbert-Schmidt (Eu- 
clidean) measure we derive tight explicit two-sided bounds for the volumes of all 
five sets. A sample consequence is the fact that, as increases, a generic positive 
map becomes not decomposable and, a fortiori, not completely positive. Due to the 
Jamiolkowski isomorphism, the results obtained for quantum maps are closely con- 
nected to similar relations between the volume of the set of quantum states and the 
volumes of its subsets (such as states with positive partial transpose or separable 
states) or supersets. Our approach depends on systematic use of duality to derive 
quantitative estimates, and on various tools of classical convexity, high-dimensional 
probability and geometry of Banach spaces, some of which are not standard. 
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1 Introduction 



Processing of quantum information takes place in physical laboratories, but it may be 
conveniently described in a finite dimensional Hilbert space. The standard set of tools 
of a quantum mechanician includes density operators which represent physical states. A 
density operator p is Hermitian, positive semi-definite and normalized. The set of den- 
sity operators of "size" 2 is equivalent, with respect to the Hilbert-Schmidt (Euclidean) 
geometry, to a three ball, usually called the Bloch ball. The set of density operators of 
"size" N forms an N"^ — 1-dimensional convex body which naturally embeds into M.Ni 
the space oi N x N (complex) matrices. 

The interesting geometry of these non-trivial, high-dimensional sets attracts a lot of 
recent attention [H El [3l HI E] . In particular one computed their Euclidean volume and 
hyper-area of their surface [6], and investigated properties of its boundary [7j. 

If the dimension N of the Hilbert space ?^ at is a composite number, the density 
operator can describe a state of a bipartite system. If such a state has the tensor product 
structure, p = Pa® Pb, then it represents uncorrelated subsystems. In general, following 
[5], a state is called separable if it can be written as a convex combination of product 
states. In the opposite case the state is called entangled and it is valuable for quantum 
information processing [9j, since it may display non-classical correlations. 

The set A4^^^ of separable states forms a convex subset of positive volume of the 
entire set of states, which we will denote by [10]. Some estimations of the relative 

size of the set of separable states were obtained in [H], [121 HSl HH UHl [lEl [I7j, while its 
geometry was analyzed in [THl [iHl [201 [21] . Similar issues for infinite-dimensional systems 
were studied in [22] . 

Quantum information processing is inevitably related with dynamical changes of the 
physical system. Transformations that are discrete in time can be described by linear 
quantum maps, or super-operators, ^ : Ai^ Aijy (or, more generally, $ : J^k 
A4n)- a map is called positive (or positivity-preserving) if any positive (semi-definite) 
operator is mapped into a positive operator. A map $ called completely positive (CP) 
if the extended map $ (g) is positive for any size k of the extension. Here Ik is the 
identity map on A^^- We will denote the cones of positive and completely positive maps 
(on A4n) by Vn and CVn respectively, or simply by V and CV if the size of the system 
is fixed or clear from the context. 

Conservation of probability in physical processes imposes the trace preserving (TP) 
property: Tr $(p) = Tr p. It is a widely accepted paradigm that any physical process may 
be described by a quantum operation: a completely positive, trace preserving map. (In 
the context of quantum communication, quantum operations are usually called quantum 
channels.) 

The set CVjf of quantum operations, which act on density operators of size N, forms 
a convex set of dimension N"^ — iV^ . Due to Jamiolkowski isomorphism [23l [2l] the set 



2 



N~^CVjf can be considered as a subset of the (A^^ — 1) -dimensional set AA^^\ of density 
operators acting on an extended Hilbert space, TIn ® 'Hn- This useful fact contributes 
to our understanding of properties the set of quantum operations, but its geometry is 
nontrivial even in the simplest case of iV = 2 p5l [25] . 

The main aim of the present work is to derive tight two-sided bounds for the Hilbert- 
Schmidt (Euclidean) volume of the set CVj^ of quantum operations acting on density 
operators of size and analogous estimates for the volume of the sets Vj^ of positive 
trace preserving maps, and of similar subsets of the superpositive cone SVn (see (1131) 
and/or [27]) or the cone of decomposable maps (see f lTTl) ) etc. We show that, for 
large N, some subsets cover only a very small fraction of its immediate superset, while in 
some other cases the gap between volumes is relatively small. These bounds are related 
to (and indeed derived from, making use of the Jamiolkowski isomorphism) analogous 
relations between the volumes of various subsets of the set of quantum states such as 
those consisting of separable states or of states with positive partial transpose (PPT) (see 
the paragraph following f|T5|) ) and their dual objects. Our methods are quite general and 
allow to produce tight two-sided estimates for many other sets of quantum states or of 
quantum maps. 

The paper is organized as follows. In the next section we introduce some necessary 
definitions involving the set of trace preserving positive maps and its relevant subsets or 
supersets, which will allow us to present an overview of the results obtained in this paper 
(summarized in Tables 2-4). Section [3] contains more definitions and various preliminary 
results. Most of those results are not new, but many of them are not well-known in the 
quantum information theory community. In section HI we state precise versions of our 
results and outline their proofs. Some details of the proofs and technical results (from 
all sections) are relegated to Appendices. 
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2 Positive and trace preserving maps: notation and 
overview of results 



2.1 Cones of maps and matrices 

Let ^ : A4n ^ -Mn Si. linear quantum map, or a super-operator. More general maps 
$ : M.K ^ M.N may also be considered and analyzed by essentially the same methods, 
but we choose to focus on the case K = N to limit proliferation of parameters. 
Let p E A4n] the transformation p' = $(p) can be described by 

where we use the usual Einstein summation convention. The pair of upper indices 
defines its "row," while the lower indices determine the "column." This agrees with 
the usual linear algebra convention of representing linear maps as matrices. The relevant 
basis of A^^r is here Eij := |ej)(ej|, i,j = l,...,N, where {ei)f^i is an orthonormal basis 
of T^AT (which can be identified with C^), and the m/i'th "column" of ^nu^, i-e., the N xN 

matrix ($ni^)^^^^, is indeed <^{E.^^,) = ^nu^Enu . 

By appropriately reshuffling elements of $rii^ we obtain another matricial represen- 
tation of a quantum map, the dynamical matrix [28], sometimes also called in the 
literature "the Choi matrix" of $. The dynamical matrix is obtained as follows 

Dmn '.= ^ nu . (2) 

mfi ^ ' 

An alternative (and useful) description of the dynamical matrix is as follows 



N 

:= {In ® $)pmax = Yl ^'"^ ® HEmf,), 

m,fi=l 



where pmax = 10(^1 5 with |^) = J2m=i ® Cm, is a maximally entangled pure state on 

We point out that the order of indices of the matrix D in ([2]) is different than in 
the previous work [2U (The reason for this change will be elucidated in the next 
paragraph.) Note that in the present notation the operation of "reshuffling," which 
converts matrix $ into D, corresponds to a "cyclic shift" of the four indices. 

It is sometimes convenient to arrange the row and column indices of D<j, and ^i, 
respectively) in the lexicographic order, thus obtaining a standard "flat" A^^ x A^^ matrix 
with a natural block structure: the leading indices m indicate the position of the block 
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and the second pair of indices n refers to the position of the entry within a block. In 
other words, the m/x'th block of D$ is or 

= (3) 

an X block matrix with each block belonging to Mn- 

If a super-operator $ belongs to the positive cone V (i.e., $ is positivity-preserving), 
then it also maps Hermitian matrices to Hermitian matrices. This in turn is equivalent to 
$ commuting with complex conjugation ^; in what follows we will generally consider only 
maps with this property. It is easy to check that Hermiticity-preserving is equivalent to 
the following relation (which has no obvious interpretation) 

^nu = ^un. (4) 

However, expressing condition (jlj) in terms of the dynamical matrix we obtain 

Dmn = D an . 

fJ-v mn 

which just means that D<p is Hermitian. Thus one may describe linear Hermiticity- 
preserving maps on via Hermitian dynamical A^^ x A^^ matrices. The property 
of being positive can be characterized just as elegantly. A theorem of Jamiolkowski [23] 
states that a map $ is positive, $ G P, if and only if the corresponding dynamical matrix 

is block positive. [A (square) block matrix (Mj^) (say, with Mij G A^at for all is 
said to be block positive iff, for every sequence of complex scalars ^ = (^j), the N x N 
matrix J2i j ^ij^i^j is positive semi-definite.] 

Arguably the most useful upshot of the dynamical matrix point of view arises in the 
study of CP maps. A theorem of Choi [29] states that a map $ is completely positive, 
$ G CV, iff Dq, is positive semi-definite. Therefore, to each CP map on A^^r corresponds 
an A^^ X A^^ (positive semi-definite) matrix, and vice versa. In particular, the rescaled 
dynamical matrix D associated with a (non-zero) CP map represents a state of a bi- 
partite system, a := D/TiD G A4*^t - see e.g. ^31 [2l], an element of the base of 
the positive semi-definite cone obtained by intersecting that cone with the hyperplane of 
trace one matrices. 

If the dimension of the cones or other sets under consideration is relevant, we will 
explicitly use a lower index, writing, e.g., CV2 for the set of one-qubit completely positive 
maps. 

2.2 Trace preserving maps 

The trace preserving property, Tr$(p) = Trp, is equivalent to a condition for the partial 
trace of the dynamical matrix 

^ Dmn = Smf, , Or TibD = I A ■ (5) 
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Therefore the compact set CVj^ of quantum operations may be defined as a common part 
of the affine plane representing the condition and the cone of positive semi-definite 
dynamical matrices - see Figure 1 in section [31 

In and (occasionally) in what follows we use the labels A, B to distinguish between 
the space on which the original state p acts, namely TIa, and the space of denoted 
1-Lb- In particular, Ia stands for the identity operator on T-La- Since such conventions are 
somewhat arbitrary (as was the ordering of indices of D), some care needs to be exercised 
when comparing ([5]) and similar formulae with other texts (such as, e.g., |26j). 

2.3 Bases of cones 

Let = {M G A^d : TrM = 0}. Next, let = {M e Md : M = c/^/^} and let 
H+ = {M e Md-^^M > 0}. If C C is a cone, we will denote by := C n H'^ the 
corresponding base of C. (This definition makes good sense if C C or, equivalently, 
if the d X d identity matrix Id belongs to the dual cone C* (see ( ITTj) ). In this case 
the cones generated by and C coincide, perhaps after passing to closures.) We will 
use the same notation for the sets of quantum maps corresponding to matrices via the 
Choi-Jamiolkowski isomorphism. Thus, for example, ^ : Ai^ ~^ -Mn belongs to 
iff TrD$ = Tr<l>(/7v) = N. (Here the identity matrix In and its image ^{In) are 
N X N matrices, while D$ is a d x d matrix, with d = N"^; in particular the two trace 
operations take place in different dimensions.) Then V fl = is a base of the cone 
V, CV n = CV^ is a base of the cone CV, and similarly for other cones that will be 
introduced later. The (real) dimension of the bases is A^^ — 1. 

The reason behind our somewhat non-standard normalization Tr M = c?^/^ is twofold. 
First, the condition can be rewritten as {M,e)Hs = !> where e = Id/d^/"^ is a matrix 
whose Hilbert-Schmidt norm is equal to one; this allows to treat e as a distinguished 
element of cones and - at the same time - of their duals. Next, the primary objects 
of our analysis are quantum maps, and the chosen normalization assures that TP (and, 
dually, unital; see Appendix 16. 5|) maps are in H^. When we are primarily interested in 
states, the normalization Tr M = 1 can be thought of as more natural (the distinguished 
element Id/d is the then the maximally mixed state, usually denoted by p^). 

While all the matrix spaces or spaces of maps are a priori complex, all cones of 
interest will live in fact in the real space M.^^ of Hermitian matrices or in the space of 
Hermi city-preserving maps. We will use the same symbols , etc. to denote the 
smaller real (vector or affine) subspaces; this should not lead to misunderstanding. 
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2.4 Other cones, all sets of interest compiled in one table 

Analogous point of view will be employed when studying other cones of quantum maps 
such as 

• the cone SP of superpositive maps (also called entanglement breaking, see f[T3|) and the 
paragraphs that follow) 

• the cone V of decomposable maps (see (fT71) ) 

• the cone T of maps which extended by identity impose positive partial transpose (see 
( lisp and the paragraphs that follow). 

In all cases we will identify the corresponding cone of A^^ x A^^ matrices and will relate in 
various ways bases of the cones and their sections corresponding to the trace preserving 
restriction. For easy reference, we list all objects of interest in the table below; see also 
Figures 1 and 3 in section [31 The missing definitions and unexplained relations (generally 
appealing to duality) will also be clarified there. 

Table 1: Sets of quantum maps and the sets of quantum states associated to them via 
the Jamiolkowski-Choi isomorphism, cf. ([6])- ([9]) below. The inclusion relation holds in 
each collumn, e.g. Vn D I^n ^ CVn D 'Tat ^ SVn- The symbols o and -k in the 
rightmost column denote sets consisting of also non-positive semi-definite matrices which 
technically are not states (and are not readily identifiable with objects appearing in the 
literature). 



Maps 






$ : Mn - 
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cones 
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-k 
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CVn 




cv% 




cvir 


KAtot 

J^In2 


PPT 
inducing 


Tn 


D 




D 




A^PPT 
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SVn 
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sv% 


D 


svir 


■/^tN2 



The action of the Jamiolkowski-Choi isomorphism, associating cones of maps to cones of 
matrices and their respective bases, can be summarized as 

$ is positive ($ G Vn) o" is block-positive (6) 
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$ is completely positive ($ G CVn) cr is positive semi-definite 

^eCV% ^ a = j^D^eMf2 (7) 

$ is PPT inducing ($ G Tyv) ^ (t G PPT 

$gT;^ ^ a = ^D^eMlV (8) 

$ is superpositive ($ G SVn) cr is separable 

<I>g5< ^ a = lz}$GA<^^P (9) 

The description of the matricial cone associated to the cone "Dat of decomposable maps 
is largely tautological: the sum of the positive semi-definite cone and its image via the 
partial transpose. We likewise have ;^-D$ G o (resp., G iff $ G (resp., G V^). 

2.5 Comparing sets via volume radii, overview of results 

Explicit formulae for volumes of high dimensional sets are often not very transparent 
(when they can be figured out at all, that is). This may be exemplified by the closed 
expression for the volume of the — 1-dimensional set A^Jf*, the set of of density 
operators of size d that has been computed in [B] 

voi(^r) = Vdi2nr^'-'y^^^^:^^. m 

Given the complexity of formulae such as (fTOl) . the following concept is sometimes 
convenient. Given an m-dimensional set K, we define vrad(i^), the volume radius of K, 
as the radius of an Euclidean ball of the same volume (and dimension) as K. Equivalently, 

vrad(i^') = (vol(i^')/vol(i?™))^''"^, where is the unit Euclidean ball. It is fairly easy 
(if tedious) to verify that (ITUl) implies a much more transparent relation vrad(A^|f*) = 
e~^^'^d~^^'^{l ± 0{d~^)) as d ^ oo, and similar two-sided estimates valid for all d. 

This point of view allows to present in a compact way the gist of our results. We 
start by listing, in Table [H bounds and asymptotics for volume radii of bases of various 
cones of maps acting on A^-level density matrices. Observe that the bounds for volume 
radii of three middle sets {T>, CV and T) do not depend on dimensionality. On the other 
hand, the volume radii of the base for the largest set V of positive maps grow as V^, 
while the volume radii of the smallest set SV of superpositive maps decrease as 1/ \/N. 

The base of the set of completely positive maps acting on density matrices of size 
is up to a rescaling by the factor equivalent to the set of mixed states A^)f* of 
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Table 2: Volume radii for the bases of mutually nested cones of positive maps which 
act on A^-level density matrices. Here rcp denotes the volume radius of the base CV^ 
of the set of completely positive maps. The last column characterizes the asymptotical 
properties, where rj^™ := limAr_»oo vrad(A'^) with X standing for V, CV or T. It is tacitly 
assumed that the limits exist, which we do not know for X ^ CV (the rigorous statements 
would involve then liminf or lim sup, cf. Theorem [S]). The question marks "?" indicate 
that we do not have asymptotic information that is more precise than the one implied by 
the bounds in the middle column. It is an interesting open problem whether r^™ admits 
a nontrivial (i.e., < 1) upper homid] cf. remark (c) following Theorem [51 



Sets of maps 


Bounds for volume radii 


Asymptotics 


positive V 
decomposable V 
completely positive CV 
PPT inducing T 
super positive SV 


\VN < vrad(P^) < 
rev < vrad(r'^) < Srcv 

\ < rcp := vrad(CP^) < 1 
Ircv < vrad(T^) < rev 
IVW < vrad(5n) < 4^ 


? 

rl"^ < 2 

„lim _ „-l/4 
'CV ~ ^ 

„lim \ 1 

rr ^ 2 
? 



dimensionality d = N'^, and similarly for other cones of maps - see eq. ([7])- ([9]). Therefore, 
the results implicit in the last three rows of Table [2] are equivalent to the following bounds, 
presented in Table [31 for the volume radii of the set of quantum states and its subsets, 
some of which were known. 

Finally, we list in Table [H the volume radii of the main objects of study in this paper: 
the set CVjf of quantum operations and of other "ensembles" of trace preserving maps. 
Each of these sets forms a N"^ — A^^ cross-section of the corresponding N"^ — 1-dimensional 
base (i.e., of CV^ etc.). 

Although the volume of the larger set is sometimes known fllOp . the cross-sections ap- 
pear much harder to analyze. Our approach does not aim at producing exact values (even 
though here and in the previous tables we made an effort to obtain "reasonable" values 
for the numerical constants appearing in the formulae). Instead, we produce two-sided 
estimates for the volume radius of CVj^ , which are quite tight in the asymptotic sense (as 
the dimension increases) and analogous bounds for the sets of positive, decomposable, 
PPT-inducing and super-positive trace preserving maps. Note that these bounds are 
similar to the results for the bases of all five sets presented in Table [21 but are not their 
formal consequences. 
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Table 3: Volume radii for the set of states A^Jf* of size d and its subsets M^^^ and A^^'^''- 
The latter two sets are well defined if the dimensionality ci is a square of an integer. Here 
a ~ 6 means that lim^^^oo a/6 = 1, while a > 6 stands for liminf(i_»oo d/h > 1. 



Sets of states 


Bounds for volume radii 


Asymptotics 


all states 
PPT states 
separable states 


< not := vrad(A^r) < 73 
iT^rtot < rpp^ := Ymd{M7^) < T^'^tot 
11 < vrad(Air) < 41 


' PPT ~ 2 
? 



While we concentrate in this work on the study of various classes of trace preserving 
maps, our approach allows deriving estimates of comparable degree of precision for other 
sets of quantum maps. As an illustration, we sketch in Appendix 16.51 an argument giving 
tight bounds for the volume of trace non-increasing (TNI) maps. An exact formula for 
that volume was recently found by a different method [30j independently from the present 
work. 

Finally, let us point out that formula ffTOj) is valid only in the case when the underlying 
Hilbert space is complex, and that our analysis focuses on the complex setting, as it is 
the one that is of immediate physical interest. However, all the discussion preceding ( fTOi) 
can be carried out also for real Hilbert spaces, and virtually all results that follow do 
have real analogues. This is because even when closed formulae are not available, the 
methods of geometric functional analysis allow to derive two-sided dimension free bounds 
on volume radii and similar parameters. Accordingly, while in the real case one may be 
unable to precisely calculate coefficients such as e~^^'^ above, it will be generally possible 
to determine the relevant quantities up to universal multiplicative constants. 
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Table 4: Asymptotic properties of volume radii for five nested sets of trace preserving 
maps. Same caveat as in Table 2 applies to the limits in the second column. Upper and 
lower bounds valid for all (as in the middle columns of Tables [2] and [3]) can be likewise 
obtained. 



Sets of trace 
preserving maps 


asymptotics of their volume radii 


positive V 
decomposable V 
completely positive CV 
PPT inducing T 
super positive SV 


\ < \im^_,^-^'p"^ < 6 
e-i/^ < limjv^oo vrad(P^P) < 2 

lim^v^oo vrad(CP^P) = e'^/^ 
i < lim^^„,vrad(rTP) < g-V^ 
1 / vr^disvin ^ A 

6 - iimTV^oo S 4 



2.6 A generic positive map acting on a high dimensional system 
is not decomposable 

This is immediate from Table HI the volume radius of the set of positive trace preserving 
maps acting on an dimensional system is of order V^, while the volume radius of the 
corresponding set of decomposable trace preserving maps is 0(1). Thus, for large A^, the 
latter set constitutes a very small part of the former one. Note that in order to compare 
volumes we need to raise the ratio of the volume radii to the power A^^ — A^^, which yields 
roughly A^~^ a fraction that is (strictly) sub exponential in the dimension of the set. 



11 



3 Known and preliminary results 



3.1 Duality of cones 

Spaces of operators or matrices are endowed with the canonical Hilbert-Schmidt inner 
product structure. The Choi-Jamiolkowski isomorphisms transfers this structure to the 
space of quantum maps. We define 

^) := {D^, D^)hs ■■= Tr dId^. 

The spaces in question and the corresponding inner products are a priori complex. How- 
ever, if we restrict our attention to the real vector spaces of Hermicity- preserving maps $ 
and Hermitian matrices which we will do in what follows, the scalar product becomes 
real and we may simply write 

We next define a duality * for cones of maps via their representation (or dynamical) 
matrix by 

C* := : Mn Mn ■■ > for all $ G C}. (11) 

This is a very special case of associating to a cone in a vector space the dual cone in the 
dual space (here Ai^ is identified with its dual via the inner product (•, ■)hs)- Duality 
for cones of matrices and cones of maps is the same by definition. 

We point out that all the cones C we consider are non-degenerate, i.e., they are of 
full dimension in the real vector space of Hermitian matrices, or in the space of 

linear maps commuting with ^ (equivalently, every map/matrix - Hermicity-prcserving 
or Hermitian, as appropriate - can be written as the difference of two elements of C) and 
further —C flC = {0}. Consequently, their duals arc also non-degenerate. 

Since the cone of positive semi-definite matrices is self-dual, it follows that 

CV* = CV. (12) 
The superpositive cone SV may be defined via duality 

Sr := V*. (13) 
By the bipolar theorem for cones ((C*)* = C), we then have 

SV* = V. (14) 
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[Note that the bipolar theorem for closed cones follows, for example, from the easily 
verifiable identity C* = —C°, where ° is the standard polar defined by K° = {x : 
{x, 2/) < 1 for all y G K}, and from the bipolar theorem for the standard polar, i.e., from 
the equality {K°)° = K valid whenever is a closed convex set containing 0.] Clearly 

SV CCV C V, 

see Figure 1. 




Figure 1: Sketch of sets of maps, a) The cone V of positive maps includes the cone CV 
of completely positive maps and its subcone CV containing the superpositive maps, dual 
to V. Trace preserving maps belong to the cross-section of the cones with an affine plane 
of dimension A^^ — A^^ (and of codimension N"^), representing the condition TtbD = J^. 
b) The sets of trace preserving maps in another perspective. This is a complete picture 
for = 2 since some of the cones coincide, namely V = V and T = SP. For > 3, the 
complete picture is more complicated, see Figure 3. 

To clarify the duality relations ([T3l) . ([T^ and the structure of the cone SV, we recall 
that $ is positive iff is block positive, which - by definition - is equivalent to > 

for every matrix of the form := that is, for every rank one positive semi-definite 

matrix. In other words, for any C, G Ti-A and for any rj ^Hb, 

0<{m){^\)v,v)HS = Trm){^\)\v){v\ = ^rD^{p^®p,) = {D^,p^^p,)Hs 

where the first tracing takes place in Hb (or J^n) and the other in Ha ® Hb, or in 
AIatz (and similarly for the two Hilbert-Schmidt scalar products). This is the same as 
saying that belongs to the cone of matrices that is dual to the separable cone (the 
cone generated by all <8) Pr? = P^^r? or, equivalently, by all products Pa® Pb of positive 
semi-definite matrices). By the bipolar theorem for cones, this is equivalent to the cone 
{D$ : $ G SV} being exactly the separable cone. 

An alternative description of SV, which justifies the "entanglement breaking" termi- 
nology, is as follows: $ is superpositive iff for every k the extended quantum map $ ® 
maps positive semi-definitive matrices to (positive semi-definite) separable matrices, or 
states to separable states if $ is trace preserving. 
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Sometimes (see, e.g., Appendix 16.21) it is useful to work with extended sets of maps 
such as the convex hulls of U {0} or V% U {0}. For technical reasons, we find the 
latter one more useful; we will denote it by = V^, and similarly for other cones. Here 
denotes the "zero" map, which may be chosen as a reference point. Further, one may 
consider symmetrized sets such as CV^^^ = CPJv™' convex hull of -CV^UCV^, where 
—CV^ is the symmetric image of CV^ with respect to 0. (Note that CV^^^ is also the 
convex hull of -CP^ U CP^, see Figure 2 below.) The advantage in using 0-symmetric 
sets is that, first, they often admit an interpretation as unit balls with respect to natural 
norms and, second, that symmetric convex bodies have been studied more extensively 
than general ones convex bodies. 







Tr(D^)=N 

















Figure 2: The set CV^ of normalized quantum maps arises as a cross-section of the 
unbounded cone of CP maps with the hyperplane representing the condition Tt D = N. 
The set CV^ of maps extended by the zero map is the convex hull of CP^ U {0}, while 
CP'^"" is the symmetrized set, the convex hull of -CV^ U CV^. CT"'^"' may be identified 
with a ball in trace class norm, whose radius equals N. Analogous notation (and similar 
identifications) may be employed for other sets of maps including V, SV etc., or for 
abstract cones. 

We next introduce the auxiliary cone of completely co-positive (CcP) maps 

CcP = {$ : To$ G CP}, 

where T : M.]^ ^ M.^ '^^ the transposition map (which is positive, but not completely 
positive for iV > 1), and the cone 

T ■.= CVnCcV. (15) 

In terms of dynamical (Choi) matrices, -Dto$ is obtained from by transposing each 
block, i.e., by the partial transpose in the second system. This means that : $ G T} 
is exactly WT, the positive partial transpose cone (positive semi-definite matrices whose 
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partial transpose is also positive semi-definite). Since, as is easy to check, separable 
matrices are in PPT, it follows that 

SV CT CCV (16) 



For N = 2 the sets T and SV coincide, while for larger dimensions the inclusion SV C T 
is proper as shown in Figure 3. 




Figure 3: Sketch of sets of maps for > 3. a) The cone V of positive maps includes 
a sequence of nested subcones: the cone V of decomposable maps, CP of completely 
positive maps, T of maps which extended by identity impose positive partial transpose, 
and the cone SP of superpositive maps, b) The sequence of nested subsets of the compact 
set of positive trace preserving maps. 

Similarly to superpositive maps, there is an alternative description of T in the lan- 
guage of extended quantum maps: $ G T iff $ (8)Xfc is PPT inducing for any size k of the 
extension, i.e., for any state p acting on the bipartite system its image, p' = $ ®I{p) G 
VVT . [The necessity of the latter condition follows by noticing that the partial transpose 
of p' equals {T ®I)p' = (T o $ (g) X)p, which is positive semi definite due to T o $ being 
CP.] 

A quantum map $ is called decomposable, if it may be expressed as a sum of a CP 
map \E'i and a another CP map \E'2 composed with the transposition T, 

$ = ^i+To^2 (17) 

or, equivalently, as a sum of a CP map and a CcP map. In other words, the cone V of 
decomposable maps is defined by 

V := CV + CcV 

(the Minkowski sum). Since the transposition preserves positivity, T) dV. It is known 
[3T|,[32] that every one-qubit positive map is decomposable, so the sets V2 and I>2 coincide. 
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However, already for = 3 there exist positive, non-decomposable maps [33], so 
forms a proper subset of V3 - see Figure 3. 

It follows from the identity ($, T o vp) = (T o $, \[/) valid for all $, \1/ that 

CcV* = CcV. 
Accordingly, the dual cone T>* verifies 

V* = {cv + Ccvy = cv* n CcV* = cvn CcV = r (is) 

This is a special case of the identity (Ci + C2)* = Q fl Q (the Minkowski sum) valid for 
any two convex cones Ci,C2- It now follows by the bipolar theorem that 

V = r*. (19) 

AsSV CT C CP by ([H]), it follows by duality that 

CP cV CP. 

3.2 Bases of cones and duality; the inradii and the outradii. 
The symmetrized sets 

We now return to the analysis of bases of cones of matrices, as defined in section 12.31 
As was to be expected, natural set-theoretic and algebraic operations on cones induce 
analogous operations on bases of cones. Sometimes this is trivial as in (CinC2)'^ = C\nC2, 
in other cases simple: (Ci + C2)^ = conv (C^UC^), where conv stands for the convex hull. 
What is more interesting and somewhat surprising is that also duality of cones carries 
over to precise duality of bases in the following sense. 

Lemma 1 Let V be a real Hilbert space, C G V a closed convex cone and let e G V be a 
unit vector such that e E CnC* . Set := {x E V : {x, e) = 1} and let and C^ = C HV^ 
and {C*y = C* nV^ be the corresponding bases of C and C* . Then 

{C*y := C*nV'' = {y eV'' : Vx G {-{y -e),x-e)<l}. (20) 

In other words, if we think of as a vector space with the origin at e, and of C^ and 
{C*y as subsets of that vector space, then (C*)*^ = —{C^y. 

Recall that for abstract cones C C V, the dual cone C* is defined (cf. ( ITT]) ) via 

C* ■.= {x eV : Vy G C {x,y) > 0}. 
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This elementary Lemma seems to be a folklore result, but does not appear in standard 
references for convexity (the best source we were pointed to after consulting specialists 
was Exercise 6, §3.4 of [M])- However, once stated, the Lemma is straightforward to 
prove. If (x, e) = (?/, e) = 1, then { — {y — e),x — e) = —{y,x) + 1 and so the condition 
from (l20l) can be restated as 

\/x eC^ - {y,x) + 1 <1 ^ \/x eC^ {y, x) > 0. 

Since under our hypotheses generates C, the latter condition is equivalent to (y, x) > 
for all X E C, i.e., to y G C*, as required. 

Let us now return to our more concrete setting of V = Ai^/ (endowed with the 
Hilbert-Schmidt scalar product) and e = Id/d}^"^- Even more specifically, we will con- 
sider V = A^^2, identified via the Choi-Jamiolkowski isomorphism with the space of 
Hermicity preserving quantum maps on TWtv; and the cones that we defined in prior sec- 
tion. Note that the quantum map associated to e = In^/N is the so-called "completely 
depolarising map," which is usually denoted by and whose action is described by 
$,(M) = (iV-iTrM)/jv. The duality relations for cones dH]), ([13]), ([14]) and ([18]), ([19]) 
combined with Lemma [1] imply now 

Corollary 2 We have the following duality relations for the bases of cones 

{cv^y = -cv^, {sv^y = -v^, (p^)° = -sv^ 

{V^y = -T^ {T^y = -V^, (21) 

where both the polarity and the negative signs refer to the vector structure in = {$ : 
Ti = Tr$(/Ar) = A^} with as the origin. 

In other words, we have for example 

V^' = {<^ eH^ : e T'^ (-($ - $,), (^ - $,)) < 1}. 

While the duality relations for cones described in the preceding subsection are rather well 
known, the duality for bases in the present generality appears to be a new observation. 
When combined with standard results from convex geometry, most notably Santalo and 
inverse Santalo inequalities [3S112S] (see below), and other tools of geometric functional 
analysis, it allows for relating volumes of bases of cones to those of the dual cones, and 
ultimately for asymptotically precise estimates of these volumes and of volumes of the 
corresponding sets of trace preserving maps. 

Let us also note here one immediate but interesting (and presumably known) conse- 
quence of the duality relations. 
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Corollary 3 For each of the sets CV^, SV^, V^, V% andT^, the Euclidean (i.e., Hilbert- 
Schmidt) in-radius is (N'^ — 1)"^/^ and the Euclidean out-radius is {N"^ — 1)^^"^. 

We observe first that, for each of the above sets, is the only element that is invariant 
under isometries of the set. Accordingly, it is enough to restrict attention to Hilbert- 
Schmidt balls centered at For CV^, the assertion is just a reflection of the elementary 
fact that A^Jf* contains a Hilbert-Schmidt ball of radius 1/ d{d — 1) centered at the 
maximally mixed state p^, and that the distance from to pure states is a/1 — 1/d. For 
SV%, it is a consequence of equality of in-radii of Ai^^l and A^^? (in the bivariate case) 
established in [13] (the out-radius of the latter is of course attained on pure separable 
states). It then follows that the in- and out-radii must be the same for the intermediate 
set T^. Finally, since the out-radius of K° is the reciprocal of the in-radius of K (and 
vice versa), we deduce the assertion for and via fl?Il) . 

It is curious to note that the statement about the out-radius of is equivalent - via 
simple geometric arguments - to the following fact (which a posteriori is true) 

IfM = {Mjk)^^j^=i is a block-positive matrix, then Tr (M^) < (TrM)^. 

It would be nice to have a simple direct proof of the above inequality, as it would yield 
(via Lemma [1] and (!2T|) ) an alternate derivation of the result from [13] concerning the 
in-radius of the set of separable states in the bivariate case. 

Similarly, the best (i.e., the smallest) constant R in the inclusion 

is the same as the best constant in 

P'^-<l>, C i?(CP'^-$,). 

It has been shown in [13] that the optimal R satisfies N"^ /2 + 1 < R < N'^ - 1. [The 
upper bound follows just from the formulae for the inradius of SV^ and the outradius 
oiCV^ (or, equivalently AI'^^p, M^"^).] A gain, there could be a more direct elementary 
argument. 

Remark 4 The Euclidean inradii and outradii of CV'Jf , SV"^ ,Vj^ ^'Djf and T^^ are 
the same as for the larger -type sets, i.e., {N'^ — 1)^^/^ and {N"^ — 1)^/^. 

As pointed out in the arguments following the statement of Corollary [31 while the fact 
that the inradii and outradii of all sets in that Corollary are identical is nontrivial, 
there is no mystery about at least some of the maps (or directions) that witness them. 
In the language of the sets of states (i.e., matrices with trace one normalization) such 
witnesses are, for outradii pure states, and universal witnesses that work for all five sets 
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are pure separable states. By duality (i.e., Lemma [T]), direction that witness inradii (for 
all sets) are obtained by reflecting a pure separable state with respect to the maximally 
mixed state p*. In the language of quantum maps purity (i.e., the Choi matrix being of 
rank one) corresponds to the map being of the form p — >■ v'^ pv (Kraus rank one), and 
the trace preserving condition is then equivalent to v being unitary. If that unitary is 
separable (i.e., a tensor product of two unitaries acting on the first and second system), 
the corresponding pure state will be separable. This means that universal witnesses of 
outradii of C'^-type sets exist also in the smaller set by the trace preserving condition 
(I5l), i.e., inside the C"'"^-type sets. Since condition ([5]) defines an affine subspace, the 
"opposite" directions giving witnesses to the inradii also belong there. 

An alternative use of duality considerations involves symmetrized sets (cf. Figure 
2). If C C V is a cone and its base, we define C^^™ := conv(— C*^ U C^); the minus 
sign referring now to the symmetric image with respect to 0. If, as earlier, e is the 
distinguished point of C fl C* defining and (C*)^, then 

(C^y"^)° = (e-C*) n (-e + C*), (22) 

where the polarity has now the standard meaning (i.e., inside the entire space V and 
with respect to the origin). In other words, the polar of C^^^ is the order interval [— e, e], 
in the sense of the order induced by the cone C* . The advantages of this approach is 
that we find ourselves in the category of centrally symmetric convex sets, which is better 
understood than that of general convex sets, and that frequently the object in question 
(C^^™ and its polar) have natural functional analysis interpretation as balls in natural 
normed spaces. One disadvantage is that in place of one very simple operation (symmetric 
image with respect to e) we have two elementary and manageable, but somewhat non- 
trivial operations (symmetrization and passing to order intervals). We postpone the 
discussion of (122]) and related issues to the Appendix. 

3.3 Volume radii and duality: Santalo and inverse Santalo in- 
equalities 

The classical Santalo inequality |35] asserts that if C is a 0-symmetric convex body 
and its polar body, then vol (K) vol (K°) < (vol(S^)) or, in other words 

vrad(K) vrad(K°) < 1. (23) 

Moreover, the inequality holds also for not-necessarily-symmetric convex sets after an 
appropriate translation, in particular if the origin is the centroid of K or of K° , a condition 
that will be satisfied for all sets we will consider in what follow. Even more interestingly, 
there is a converse inequality [36], usually called "the inverse Santalo inequality," 

vrad(K) vrad(K°) > c (24) 
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for some universal numerical constant c > 0, independent of the convex body K (sym- 
metric or not) and, most notably, of its dimension m. 

The inequalities (125]) . fl24l) together imply that, under some natural hypotheses (which 
are verified in most of cases of interest), the volume radii of a convex body and of its polar 
are approximately (i.e., up to a multiplicative universal numerical constant) reciprocal. 
By Lemma [H the same is true for the base of a cone and that of the dual cone. This 
observation reduces, roughly by a factor of 2, the amount of work needed to determine 
the asymptotic behavior of volume radii of, say, sets from the third column of Table 1. 
We note, however, that since, at present, there are no good estimates for the constant 
c from ^QM if K is not symmetric, it is often more efficient to revisit arguments from 
[16] which allow to estimate volume radii of polar bodies without resorting to the 
inverse Santalo inequality. (An argument yielding reasonable value of c for symmetric 
bodies was recently given in [TF].) 
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4 Volume estimates: precise statements and approx- 
imate arguments 

The results stated in section [3] allow us, in combination with known facts, to determine the 
asymptotic orders (as N — >■ cx)) for the volume radii (and hence reasonable estimates for 
the volumes) of bases of all cones of quantum maps discussed up to this point. Our goal 
is slightly more ambitious; we want to find not just the asymptotic order of each quantity, 
but also establish inequalities valid in every fixed dimension and involving explicit fairly 
sharp numerical constants. Specifically, we will show the following 

Theorem 5 We have the following inequalities, valid for all N, and the following asymp- 
totic relations 

(i) i < vrad(CP^) < 1, lim^^^o vrad(CP^) = e"^/^ ^ 0.779 

(ii) iN^^^ < vrad(P^) < 6N^/^ 

(iii) < vrad(5P^) < AN~^/^ 

f \ 1 / "ad(rb) vrad(rb) 

(iv) j< / \ < 1, ^<limmfAr / \ 

V 1 < < 8, limsup,v < 26^/^ 

vrad(CP^j vrad (CP^j 

Remarks : (a) Estimates on volume radii listed in Table [2] are either identical to the 
corresponding inequalities stated above, or follow by the same argument. 

(b) Since the asymptotic orders of the volume radii of the families CV%, and are 
the same, we chose - for greater transparence - to compare the volume radii of the two 
latter sets to that of CV^ in (iv) and (v), rather than give separate estimates for each of 
these quantities. 

(c) It is an interesting open problem whether there exists a universal constant a < 1 
such that vrad(r^) < avrad(CP^) for all > 2 or, equivalently, "is vol(A^^^^) < 
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vo\{M'°l) for some a < 1 and all N > 2?" Analogous question may be asked about 
comparing vrad(P^) and vrad(C'P^). Inquiries to similar effect can be found in the 
literature [IDIES]. 

(d) It is likely that the asymptotic bound 2e^/^ in (v) holds actually for all A^. Indeed, 
there is a strong numerical evidence that the estimate vrad(C'P^) > e~^/^ from (i) is 
valid for all A^ and not just in the limit. (In view of the explicit character of the formula 
([To]) this issue shouldn't be too difficult to resolve.) Should that be the case, the next step 
would be to carefully analyze the dependence of vrad(r'^) on A^ given by the arguments 
presented in this paper. 

Since the bases of cones, whose volume radii are described by Theorem 0, are effec- 
tively homothetic images, with ratio A^, of the corresponding sets of trace one matrices 
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(see Table 1 and the formulae that follow it), some of the inequalities/relations of Theo- 
remOfoUow from known estimates for the volumes of various sets of states, particularly if 
we do not insist on obtaining "good" numerical constants that are included in the state- 
ments. For example, the estimates in statement (iii) are contained in Theorem 1 from 
[16]; one obtains the constants | and 4 by going over the proof of that Theorem specified 
to bilateral systems. Similarly, the statement (iv) is (essentially) a version of Theorem 
4 from [H] which asserts that, in the present language, vrad(A^^^"'")/vrad(A1^2) > cq 
for some constant Cq > independent of the dimension N (the upper estimate with con- 
stant 1 is trivial). However, the argument from [16] yields only Cq = | and for the 
asymptotic lower bound. 

Next, the asymptotic relation in (i) follows from the explicit formula (ITOl) : see the 
comments following (ITOl) . Presumably, the estimates in (i) can also be derived from (JTOl) . 
but there are more elementary arguments. For a simple derivation of the lower bound from 
the classical Rogers- Shephard inequality [39] see [1^, section II. And here is an apparently 
new proof of the upper bound: combine the duality results of the preceding section, 
specifically the identification {CV^)° = —CV^ from (12T1) . with the Santalo inequality (l23i) 
to obtain 

1 > vrad(CP^)vrad((CP^)°) = vrad(CP^)vrad( - CV^) = vrad(CP^)^ 

as required. We recall that, in the context of (I2T]) . the operations ° and — take place in 
the space of quantum maps verifying TrD$ = A^, with thought of as the origin; 
note that is the centroid of CV^ and so (l23l) with K = CV^ indeed does apply in that 
setting. 

Arguments parallel to the last one lead to versions of the remaining statements with 
some universal constants. For example, the identification {SV^)° = —V^ combined with 
the Santalo inequality ( l23l) and its inverse ( 12^ leads to 

1 > vrad(5P'')vrad(P^) > c, 

where c is the (universal) constant from (!24|) . Combining the above inequality with 
(iii) we obtain ^N^/^ < vrad(P^) < QN^'"^. Similarly, 1 > vrad(T^)vrad(r'^) > c 

vrad 1 I 

combined with (the already shown version of) (iv) and with (i) implies ) < 32 

vrad \CV\\ 

vrad f I 

and limsup/v } ^ \ < 4e^'^. As the constants in Theorem [5] are not meant to be 

vrad (CP^j 

optimal, we relegate the somewhat more involved (but still based on classical facts) 
arguments yielding them to Appendix 16.11 

The inequalities of Theorem [S] compare volumes of bases of cones, that is, sets of 
maps $ normalized by the condition that the trace of the corresponding Choi (or 
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dynamical) matrix, is N (or Tr$(Jiv) = A^). [Of course, any other normalization - 
most notably Tr = 1 leading to sets of states - would work just as well for comparing 
volumes provided we were consistent.] However, if we want to study quantum operations, 
i.e., trace-preserving quantum maps (or, similarly, unital maps), then - as explained in the 
previous sections - the corresponding constraints are stronger than just normalization by 
trace: in each case we are looking at an A^^-codimensional section of the cone as opposed 
to the 1-codimensional base. However, in either case the codimension is much smaller 
than the dimension, which is A^^ — N"^. The following technical result will imply that 
then, under relatively mild additional assumptions assuring that the base of the cone is 
reasonably balanced (which will be the case for all the cones we studied), the volume 
radius of the section will be very close to that of the entire base. 

Proposition 6 Let K be a convex body in an m- dimensional Euclidean space with cen- 
troid at a, and let H be a k-dimensional affine subspace passing through a. Let r = tk 
and R = Rk be the in-radius and out-radius of K. Then 

/ \ — I 

(vrad(ir) R~^ b{m, A;) V < vrad(K n H) < i VTad{K) b{m, k) 

(25) 

where b[m, k) := — , , . ^ — , ' ..^ 




The proof of the Proposition is relegated to Appendix 16. 3j now we explain its conse- 
quences. First, let us analyze the parameters that appear in fl25l) . By Corallary [3l 
for all bases of cones that we consider here we have r = — 1 and R = a/c/ — 1, 

where d = N"^. Next, we have m = d"^ - 1 = - 1, k = d"^ - d = - N"^ and 
m-k = d- l = N'^-l, in particular f = l + ^ = l + ^ and ^ = = 

Further, the quantity 6(m, k) = j " (related to the Beta function) 

is easily shown to satisfy 1/v^ < b(m,k) < 1 (for our values of m,k it is actually 
1 - C'(M)). Similarly, 1 < (™)" < 2 for all k,m and 1 + 0(M) for our values of 
m, k. Consequently, if vrad(i^') is subexponential in d (in our applications it is a low 
power of A^, hence of d), then yYa.d{K fl if)/vrad(i^r) — 1 as A^ — >• oo. 
This leads to 



Theorem 7 We have the following asymptotic relations 
(i) limjv^„oVrad(CP^P) = e'^^ 

(n) J < hmmfAT — < hmsup^ — jyi/^ ^ ^ 
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7. < hmmfjv — ^fVi/a < limsup^ — J- 1/2 < 4 



^_i/2 ^ 11111 oup TV Ar-1/2 

iv) ^ < lim mf AT V " / 



vrad 



(v) limsup^^^^4^<2e^/' 

^ ^ vrad(cP5^) - 

Upper and lower bounds in the spirit of Theorem [5] (i.e., vahd for all A^) can be likewise 
obtained. 

The reader may wonder why we perform our initial analysis on bases of cones rather 
than working directly with the smaller sets of trace preserving maps. The reason for this 
is two-fold. First, the bases being homothetic to various sets of states, any information 
about them is at the same time more readily available and interesting by itself. Second, 
while we do have - as a consequence of Lemma [T] - nice duality relations between bases 
of cones, similar results for sets of trace preserving maps are just not true. As a demon- 
stration of that phenomenon we show in Appendix 16.41 that . in contrast to the bases CV^, 
the sets CV^^ are very far from being self-dual in the sense of fl2Tl) . 



5 Conclusions 

We derived tight explicit bounds for the effective radius (in the sense of Hilbert-Schmidt 
volume), or volume radius, of the set of quantum operations acting on density matrices 
of size A^, and for other convex sets of trace preserving maps acting such matrices such 
as positive, decomposable, PPT inducing or superpositive maps. The novelty of our 
approach depends on systematic use of duality to derive quantitative estimates, and on 
technical tools, some of which are not very familiar even in convex analysis. 

Since the volume radii of the sets of trace preserving maps that are positive display a 
different dependendce on the dimensionality than those of the smaller set of decomposable 
maps, the ratio of the volumes of the latter and the former set tends rapidly to as the 
dimension increases. In other words, a generic positive trace preserving map is not 
decomposable and, a fortiori, not completely positive. Thus we were able to prove a 
stronger statement than the one advertised in the title of the paper. Similarly, a generic 
PPT inducing quantum operation (and, a fortiori, a generic quantum operation) is not 
superpositive. Analogous relations (some of which were known) exist between the sets of 
states related to those of maps via the Jamiolkowski isomorphism. 
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6 Appendices 

6.1 Better constants in Theorem [5t mean width, Urysohn in- 
equaUty and related tools 

The arguments given in the preceding section did not yield the asserted values of the 
constant | in part (ii), the constants and j in part (iv), and the constants 8 and 
in part (v) of Theorem O We will now present the somewhat more involved line of 
reasoning that does yield these constants. 

The following concepts will be helpful in our analysis. If C is a convex body 
containing the origin in its interior, one defines the gauge of K via 

\\x\\k := inf{t > : x e tK}. 

Roughly, \\x\\k is the norm, for which K is the unit ball, except that there is no symmetry 
requirement. Next, the mean width of K (or, more precisely, the mean half-width) is 
defined by 

w{K) := / ||x||i<:o(ix = / max{x, y)dx 

(integration with respect to the normalized Lebesgue measure on S"^'^). A classical 
result known as Urysohn's inequality (see, e.g., [H]) asserts then that 

vrad(K) < w{K). (26) 

A companion inequality, which is even more elementary, is 

YTad{K) > w{K°)-\ (27) 

The proof of (1271) is based on expressing the volume as an integral in polar coordinates and 
then using twice Holder inequality: vrad(_ft') = {Jgm^i W^WK^dxY^"^ > Jgm-i \\x\\]^^dx > 
(/gm-i ||a;||i^(ia;) ^ = w{K°)-^. 

Applying fl26l) in our setting of the A^^ — 1-dimensional space and for K = T>\, 
we obtain 

vrad(P^) < w{V\) = w{conY{CV^^UCcP^)) < wiCV^ + CcV^) = 2w{CV\) < 4, (28) 

because w{-) commutes with the Minkowski addition (of sets), and because w{CV%) = 
w{CcP^) < 2. The latter is a consequence of similar estimates for the set of all states 
(which is equivalent to CV^ up to a homothety), see [TH (TB]. (We note that while the 
limit relation w{CV\) ^ 2 as A^ ^ cxo follows easily from well-known facts about random 
matrices, the estimate valid for all N requires finer arguments such as those presented in 
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appendices of [H]). Combining the above estimate with part (i) of Theorem 3 we obtain 
the upper estimate in part (v) with the asserted constant 8. 

The same bound w(V^) < 4 combined with (127|) (apphed this time with K = T^) 
and with part (i) leads to the lower bound j in part (iv). 

To obtain the asymptotic bounds from part (iv) and (v) with the required constants 
^ and 2e^/^ we argue similarly, but instead of the universal estimate w{T>^) < 4 we 
use a tighter asymptotic bound limsup^ i(7(P^) < 2 This bound is a consequence of 
classical isoperimetric inequalities and the measure concentration phenomenon that they 
induce (see, e.g., [IQ]): a Lipschitz function on S"^~^ is strongly concentrated around its 
mean. In particular, if the out-radius of K is at most R, then Jgm-i \ ll^^H^ft-o — w{K) \ dx = 
0{R/m^l'^). \{ K = CV% or CcP^, then, by Corollary El R = {N'^ - l)^/^ ^j^jjg ^ ^ 
A^^ — 1, hence R/m}/'^ = (N'^ + 1)~^/^ < A^^-*^. It is then an elementary exercise to show 
that 




< max {w{CV^^),w{CcV%)) + 0{N-^) 

< 2 + 0{N~'), 

whence limsup^ 'w;(P^) < 2, as required. Universal (as opposed to asymptotic) upper 
bounds on vrad(P^) better than 4 obtained in (125]) can also be derived this way, most 
efficiently by converting spherical integrals to Gaussian integrals and using the Gaussian 
isoperimetric inequality. (This would also improve somewhat the bounds ^ and 8 in 
parts (iv) and (v), but we will not pursue this direction here as the payoff doesn't seem 
to justify the effort.) 

Finally, to obtain the lower bound on vrad('P^) from part (ii) of the Theorem, we 
note that the upper bound 4A^~^/^ for vrad(iS'P^) (stated in part (iii)) was de facto (see 
[16j) deduced from the stronger estimate w{SV%) < AN~^^'^. It then remains to apply 
(1271) and the duality between and SV%. 

6.2 Symmetrized bodies and order intervals 

We will now analyze the polar of the symmetrized body C^^™. Recall the notation of 
section 3.2: V is a real Hilbert space, C C V a closed convex cone, C* the dual cone. 
Next, e G C n C* is a unit vector, := {x E V : {x, e) = 1} is an affine subspace of V 
and = C nV^ is the base of the cone C. Finally, the symmetrized body is defined as 
Qsjm ._ conv(— UC^); the minus sign referring to the symmetric image with respect to 
0. An important point, following from classical results [33111] and explained in Appendix 
C of p!lj, is that under mild assumptions which are satisfied for all the cones we consider, 
the volume radii of and of C^^™ differ by a factor smaller than 2. 
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Our main assertion (equation (122!) in section 3.2) is that 

(C«y-)° = (e-C*) n i-e + C*), 

where the polarity has the standard meaning (i.e., inside the entire space V and with 
respect to the origin). That is, y G (C**^™)° iff both y + e and e — y are in C* or, in other 
words, iff y belongs to [— e, e], the order interval in the sense of the order induced by the 
cone C*. For example, if we want to investigate and P'^y™ = conv (— U P^), we 
may specify the framework above to C = V, obtaining 

(psym^o ^ _ _^ p (^^^ _ _ 

To prove the assertion, denote V~ := {x G V : {x,e) < 1} (one of the half-spaces 
determined by K^) and = CnV- (cf. Fi gure 2 in section [3]). Then 

C'^y"! = conv(-C^ U C^) = conv (-C^ U C^) . 

Hence, using standard rules for polar operations (see, e.g., IH]), 

(C^n°= (-c^)°n (c^)°. 

Next, 

{C^y = {CnV-y = comr{{V-yu C°) = cohv((-oo, l] ■ e U -C*) = e-C*, 

where the bar stands for the closure. Combining this with the preceding formula and 
again using the standard rules gives 

(C^y-)° = (e-r)n(-e + r) 

or the intersection of two cones with vertices at e and — e. Clearly this does not equal 
^Q*yym gxcept in dimension 1. However, the two bodies are closely related. For example, 
if e is the point of symmetry of C^, then (C*)*'^™ is a cylinder with the base {C*y and 
the axis [— e,e], while (C^y™)" is a union of two cones whose common base is (C*)^ — e, 
the central section of the cylinder, and the vertices are — e and e. The two bodies only 
differ in one dimension; if thought of as unit balls with respect to the corresponding 
norms, the two norms coincide on the hyperspace Vo := {x G V : (x, e) = 0} and 
on the complementary one-dimensional space Me, but on the entire space we have in 
the first case the direct sum in the i^o sense, while in the second case in the ii sense. 
If the base is non-symmetric, the situation is more complicated. For example, the 
section Vo H (C^™)° is congruent to the intersection of (C*)^ with its symmetric image 
with respect to e, but (see |42j) the volume radii of the two bodies are comparable if, for 
example, e is the only point that is fixed under isometrics of {C*y (as is the case in all 
our apphcations), or just the centroid of (C*)^. 
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6.3 Proof of Proposition 2: for "balanced" cones, and 
have comparable volume radius 

We may assume that a = (otherwise consider K — a). By hypothesis, we have then 

r5™ dK C. RB^, (29) 

where is the m-dimensional unit Euchdean ball. For a subspace E, denote by Pe the 
orthogonal projection onto E. Then (see |l2l H3]). 

volmiK) <YohiKnH) yoKPh^K), (30) 

where s = m — k and is the m — fc-dimensional space orthogonal to the /c-dimensional 
subspace H. Therefore 

voU(K) ^ Yoh{KnH) yoI{Ph±K) volfc(5|)vol,(5|) 



YolUBT) ~ V0lfc(5|) VOl,(5|) Y0\m{Bir) 

Hence, using f l29l) . 

vrad(ir)- < vrad(irni/)^i?^"°^^^^^')"°^^^^^^^ 



V0l„(52"^) ' 

which is the first inequality in fl25|) . For the second inequality, we start with the even 
more classical result (see [H] or [15]; same notation as fl30|) ) 

m 



VoU(ir) > (^^J Y0\k{KnH)Y0\s{PH^K), (31) 

which doesn't even require that H passes through the centroid of K. As above, this can 
be rewritten in terms of volume radii as 

7Vrad(Kr > vrad(irng)^ voU(i?p vol (i^l) 

kj ^ ^ - ^ ^ V0l,(52™) 

which is the second inequality in (J25i) . 



6.4 "No duality" for CVY 

The purpose of this Appendix is to show that, in contrast to the bases CV^, the sets CV"^^ 
are very far from being self-dual in the sense of fl?Il) . that is, that the polar of CV^^ inside 
the space defined by the trace preserving condition (jS]) considered as a vector space with 
as the origin is quite different from the reflection of CV^^ with respect to 



29 



Generally, if K C R*" is a convex body containing the origin in its interior and 
H C MJ^ is a vector subspace, K° fl if is always contained in the polar of K (1 H inside 
H, and the discrepancy between the two (i.e., the smallest constant A > 1 such that 
the polar of fl if is contained in X{K° n H)) is the same as the discrepancy between 
K r\H and the orthogonal projection of K onto H. That discrepancy is also equal to the 
maximal ratio between 

max(-u, x) and max (-u, y) (32) 

x&K y£KnH 

over nonzero vectors u E H . 

In our case K = CV% and K (1 H = CVj^. As a vector space, H may be identified 
with maps whose dynamical matrix has partial trace equal to 0. We will argue in the 
language of dynamical (Choi) matrices considered as "fiat" block matrices. In these 
terms, membership in H is equivalent to each block being of trace 0. We will choose as u 
the block matrix whose 11-th block is U = En — N~^I]sf and the remaining blocks are 0. 
Further, we will choose as x the matrix whose 11-th block is X = NEu and the remaining 
blocks are 0; then the scalar product corresponding to {u, x) is ti{UX) = N — 1. On the 
other hand, if Y is the 11-th block of the Choi matrix of any element of CVj^, then y is a 
state and so the scalar product corresponding to {u, y) is tiiUY) = ti^EuY) —N~HyY < 
tiY — N~^tTY = 1 — A^~^. Accordingly, the discrepancy between the two maxima in fl32|) 
is at least {N - 1)/(1 - N'^) = N. 



6.5 Volume radius of the set of trace non increasing maps 

We want to determine the asymptotic order of the volume radius of the set of all com- 
pletely postive, trace non increasing maps ^ : M.n ^ -Mn i-e. the set 

As pointed out earlier, an exact formula for that volume was very recently found (in- 
dependently from this work and by a different method) in [30]. However, an argument 
using the approach of this paper is conceptually very simple and so we include it. We 
have 

Proposition 8 We have, for all N, 

(e AT^/^) < vol (CP™) ^ .33^ 

^ ^ ~ Yo\{CV'Jf) X vo\{{M e Mn ■■ < M < In}) ~ 

To derive estimates on vol(C'P™^) from the Proposition, one needs to use the readily 
available information on the two factors in the denominator of the middle term of ( 133|) . 
First, the asymptotic order of the volume radius of CVj^ was determined in Theorem 
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[7(i). Next, the set A := {M E Mn ■■ < M < In} is a ball of radius 1/2 (in the 
operator norm) centered at In^^, and so its volume radius admits easy bounds given by 
the in- and outradius: 1/2 and \/N /2 (actually a much tighter lower bound \/iV /4 can 
be obtained via a slight modification of the argument from Theorem 0(1) , see Appendix 
16.11 but for our purposes the trivial bounds suffice). A straighforward calculation leads 
then to 

Corollary 9 

hm vrad(CP™) = e'^/^ 

The key point is that in order to calculate the volume radius we need to raise the volume 
to the power Thus the factors such as (e N^^^)~^^ on the left hand side of (1551) 

are inconsequential since it leads to an expression of the form 1 - 0(M). por the 
same reason, the effects of yol{A) and of the b{m, fc)-type factor, which also enters the 
calculation (cf. Proposition [6] and the comments following it), tend to as ^ oo. 

For the proof of Proposition [8] we note first that CV™^ is canonically isometric to the 
set of suhunital maps 

CV"";^ ■.= {^eCV: Wn)<In}. 

In what follows we will work with the latter set. The isometry, which assigns to $ : 
M.N ^ M-N the dual (in the hnear algebra, or Banach space sense) map $*, sends CV^ 
to the set of unital maps CV^ := {$ G CVn '■ ^(In) = In}- The set CV^ admits a 
natural fibration: with every M G we may associate 

Fm = {$ G CP : <^{In) = M}; (34) 

in particular Fj^^ = CV^. In the language of Choi (dynamical) matrices the equality from 
flM|) translates to ti = M, or to Y2j Djj = M if we think of D<^ as a block matrix 
Di^ = (Djk) = [^(Ejk)). Since all fibers Fm are parallel to the subspace M defined by 
ti bD,}, = (or J2j Djj = 0)) can express the volume as an integral 

vol(CP™) = vo\{CVf) = Ar-^V2 I ^35) 

J A 

The reason for the factor iV~^^/^ is that while the fibration is naturally parametrized by 
the elements of A, the projection of Fm onto A/""*" is actually the map p — > A^~^tr(p)M, 
whose Choi matrix is N~^In (8> M (or a block matrix whose all diagonal blocks are M/N 
and off-diagonal blocks are 0). Now, the Hilbert-Schmidt norm of N'^I^ is A^~^/^, and 
so the projection of CV^ onto A/""*" is isometric to N'^^'^A. 

The second inequality in fl551) is now an immediate consequence of ([35]) and the 
bound vol (Fa/) < vo1(F/^) = yo\{CV^) = vo\{CVjf), vahd for all M e A, which 
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in turn follows, e.g., from Fm being the image of F/^^ under the contraction qm ■ $(•) — > 
M2$(-)M2. (On the level of dynamical matrices, the action of qm is given by 
{In®M'2^D^{Im®M2^ or, in the language of block matrices, by (-Djfc) — * (M^DjkM^).) 
The fact that gM{Fij^) C Fm is obvious from the definition; surjectivity for invertible 
M's follows by considering the inverse = Qm-^i and for singular M's by looking at 
invertible approximants. 

For the first inequality in (l33l) we may use (|3T|) with K = CV^^ and the section K n 
H = CV^. As pointed out earlier, the set PhK = Pj^± is then isometric to N'^^'^A, and 
it remains to use the elementary bound ('^) = (™) < (em/s)^. A alternative argument 
is to restrict the integration in fl35|) to {M : tl^ < M < In}, which is a ball in the 
operator norm of radius (1 — 1)/2, then use the fact that for such M the function cf/^/ is 
a contraction, and finally optimize over t G (0, 1). This approach allows in fact to express 
the Jacobian of qm in terms of eigenvalues of M and, subsequently, to express the ratio 
under consideration as a multiple integral over [0, 1]^, but we will not pursue this path 
further. 
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